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ABSTRACT 

We  compare  the  performance  of  three  heuris¬ 
tic  algorithms  [3,6,13]  for  the  minimization  of 
sum-of-products  expressions  realized  by  the  newly 
developed  multiple-valued  programmable  logic 
arrays  [9].  Heuristic  methods  are  important 
because  exact  minimization  is  extremely  time- 
consuming.  We  compare  the  heuristics  to  the 
exact  solution,  showing  that  heuristic  methods 
are  reasonably  close  to  minimal.  We  use  as  a 
basis  of  comparison  the  average  number  of  pro¬ 
duct  terms  over  a  set  of  randomly  generated 
functions.  All  three  heuristics  produce  nearly  the 
same  average  number  of  product  terms. 
Although  the  averages  are  close,  there  is  surpris¬ 
ingly  little  overlap  among  the  set  of  functions 
where  the  best  realization  is  achieved.  Thus, 
there  is  a  benefit  to  applying  different  heuristics 
and  then  choosing  the  best  realization. 


I.  INTRODUCTION 

The  minimization  of  sum-of-products  expres¬ 
sions  in  binary  logic  has  received  considerable 
attention  for  over  30  years.  The  complexity  of 
the  problem  has  been  known  for  almost  as  long. 
Although  minimal  sum-of-products  extraction 
appears  to  be  a  special  case  of  the  general 
minimal  set  covering  problem,  it  is  not.  This  was 
proven  by  Gimpel  [7]  in  1965.  That  is,  any 
instance  of  the  set  covering  problem  is  an 
instance  of  the  minimal  sum-of-products  extrac¬ 
tion.  In  1972,  Karp[8]  showed  that  the  set  cover¬ 
ing  problem  is  NP  complete;  thus,  so  also  is 
minimal  sum-of-products  extraction.  The  best 
known  algorithm  then  requires  exponential  time. 


This  is  a  real  barrier;  it  precludes  the  exact 
minimization  of  functions  with  even  a  moderately 
low  number  of  inputs,  e.g.  20.  As  a  result,  con¬ 
siderable  effort  has  been  devoted  to  heuristic 
minimization  methods.  For  example,  among  the 
Berkeley  VLSI  tools  is  ESPRESSO-IIC  [4],  a  C 
program  that  minimizes  binary  functions  by  a  set 
of  operations  on  the  prime  implicants. 

Recently,  there  has  been  considerable  interest 
in  multiple-valued  PLA’s  [1,2,3,5,6,9,10,11,13,14, 
15).  Several  have  been  proposed  [11,14,15]  and  at 
least  one  implemented  [9].  We  know  of  three 
heuristic  MVL  sum-of-products  minimization 
algorithms.  Pomper  and  Armstrong  [13]  intro¬ 
duced  in  1981  a  heuristic  method  that  found  a 
near-minimal  sum-of-products  expression  as  a 
direct  cover  of  the  function.  The  algorithm 
proceeds  in  two  steps;  1.  select  a  minterm  and  2. 
find  an  implicant  that  covers  the  minterm.  The 
first  step  is  accomplished  by  a  random  choice  of 
minterm,  while  the  second  step  is  accomplished 
by  choosing  the  largest  implicant.  In  1986, 
Besslich  [3]  introduced  another  direct  cover 
method  that  seeks  to  cover  the  "most  isolated" 
minterms  first.  And  in  1987,  Dueck  and  Miller 
[6]  introduced  a  method  which  also  seeks  the 
most  isolated  minterm  first,  but  chooses  a  pro¬ 
duct  term  that  tends  to  introduce  the  fewest 
discontinuities  when  subtracted  from  the  func¬ 
tion. 

There  has  been  little  study  of  the  relative 
merits  of  heuristic  algorithms.  To  the  credit  of 
Brayton,  Hachtel,  McMullen  and  Sangiovanni- 
Vincentelli  [4],  the  realizations  produced  by 
ESPRESSO-IIC  were  compared  with  the  realiza¬ 
tions  of  MINI,  PRESTO,  and  POP  [4]  over  a  set 
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of  56  specifically  chosen  binary  functions.  How¬ 
ever,  the  justification  of  the  newly  introduced 
algorithms  examined  by  us  has  rested  on  an 
intuitive  notion,  supported  by  examples.  In  this 
paper,  we  analyze  six  synthesis  methods. 

1.  Random  minterm  random  implicant 

2.  Pomper  and  Armstrong  [13] 

3.  Besslich  [3] 

4.  Dueck  and  Miller  [6] 

5.  Gold 

6.  Absolute  minimization 

over  a  set  of  7000  randomly  chosen  2-variable  4- 
valued  functions.  Our  choice  of  2-variable 
functions  is  influenced  by  a  need  to  compare 
heuristics  (1-5)  with  absolute  minimal  realizations 
(6).  Even  with  only  2  inputs,  absolute  minimiza¬ 
tion  requires  high  computation  times,  on  the 
order  of  days.  Our  choice  of  4-values  is 
influenced  by  the  current  interest  in  this  radix. 
The  three  heuristic  algorithms  in  [3,6,13]  apply  to 
sum-of-products  expressions  where  sum  is  MAX. 
In  all  three  algorithms,  the  case  where  sum  is 
truncated  SUM  is  also  considered.  Because  the 
PLA’s  produced  in  CCD  [9]  use  the  SUM  opera¬ 
tion,  our  analysis  uses  this  operation.  Further, 
we  have  adapted  the  heuristics  to  the  case  of  pro¬ 
duct  terms  consisting  of  interval  literals,  again 
since  the  CCD  PLA’s  realize  such  functions.  A 
motivation  for  the  study  reported  here  is  the 
development  of  a  minimization  method  for  our 
CCD  PLA  CAD  tool  [10]. 

This  paper  is  organized  as  follows.  The  next 
section  introduces  notation  and  fundamental  con¬ 
cepts.  Section  HI  describes  the  six  synthesis 
methods,  and  Section  IV  describes  the  results  of  a 
comparative  analysis  of  their  performance.  In  the 
final  section,  we  summarize  the  results. 


II.  BACKGROUND  AND  NOTATION 

Let  X  =  {a:j,  *2  >  ‘  ‘  ^  ”  vari¬ 

ables,  where  a:,-  takes  on  values  from 
R  =  {0, 1, ...,  r— 1}.  A  function  f  {X)  is  a  map¬ 
ping  f:R"—*R  lj{r},  where  r  is  the  don’t  care 
value.  Specifically,  /  (X)  is  said  to  be  an  n- 
variable  r-valued  function.  Fig.  1  below  shows  a 
map  representation  of  a  2-variable  5-valued  func¬ 
tion.  Blank  entries  correspond  to  0.  An  assign¬ 
ment  X  of  values  to  variables  in  X  is  called  a 
minterm  if  /  (a;)  0.  In  Fig.  1,  there  are  12  min- 

terms  all  of  which  yield  a  2. 


Functions  realized  by  the  PLA’s  described  in 
Kerkhoff  and  Butler  [9]  are  composed  of  three 
functions, 

1.  MIN:  f{xi,X2)  =  X1X2  (=  MIN[xi,xf)), 

2.  SUM:  f{xi,X2)  =  X1+X2  (=  x^+X2  if  X1+X2 
<  3  and  =  3  otherwise,  where  -I-  is  viewed  as 
ordinary  addition  and  a:,-  is  viewed  as  an 
integer),  and 

3.  literal:  /(xj)  =  “xj  (=  r— 1  if  a  <Xj<6 
and  =  0,  otherwise). 


2  411 


2  X 


1  *2 


Figure  1.  Example  of  a  2-variable  5-valued  function. 


In  binary,  the  SUM,  MIN,  and  literal  correspond 
to  AND,  OR,  and  x* ,  where  x*  €  {x,x}.  In  the 
realization  of  functions  by  a  multiple-valued 
PLA,  constants  and  literals  occur  as  operands  of 
the  MIN  functions.  A  product  term  is  the  MIN 
of  one  nonzero  constant  and  one  or  more  literals. 
For  example,  f  {xi,x,f)  =2^Xi  ^X2  is  a  product 
term  that  is  2  when  Xj  is  2,  3,  or  4  and  X2  is  1. 
An  implicant  of  function  /  (X)  is  a  product  term 
/(X)  such  that  /(X)>/(X).  A  prime  implicant 
of  /(X)  is  an  implicant  of  /(X)  such  that  there 
is  no  other  implicant  t  (X)  of  /(X)  where 
/'(X)>/(X).  For  example,  2^Xj  X2  is  an 
implicant  of  the  function  in  Fig.  1.  but  it  is  not  a 
prime  implicant.  However,  2  ^x^^  ^X2  is  a  prime 
implicant.  Any  function  can  be  expressed  as  the 
SUM  of  implicants  [15].  For  example,  the  func¬ 
tion  f{xi,X2)  in  Fig.  1  can  be  expressed  as  the 
SUM  of  4  implicants. 


2^x/  ^Xj  +2^xf  “^x^  +2^x^  ^x^  +2^x}  ^x^.  (1) 


We  use  the  term  sum-of-products  to  describe 
functions  realized  by  such  PLA’s,  where  sum 
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refers  to  SUM  and  product  refers  to  MIN.  A 
sum-of-products  expression  for  function  f  {X)  is 
minimal  if  there  is  no  other  expression  for  /  (X) 
with  fewer  product  terms.  For  example,  (1)  is  a 
minimal  sum-of-products  expression.  Given 
/  (X),  implicant  /(X)  covers  a  minterm  at  x  if 
/  lx)  =  lix).  Therefore,  !7(X)  =  /  (X)  -  /(X) 
has  the  property  g{x)  =  0,  and  we  say  that  sub¬ 
tracting  /(X)  drives  the  minterm  at  x  to  0. 

It  is  interesting  to  note  that  the  none  of  the 
implicants  in  (1)  are  prime.  In  fact,  any  sum-of- 
products  expression  for  /(a:i,a:2)  containing  at 
least  one  prime  implicant  is  not  minimal.  Such 
expressions  require  at  least  5  product  terms. 
Therefore,  in  contrast  to  the  case  of  conventional 
binary  sum-of-products,  it  is  necessary  to  con¬ 
sider  all  implicants.  A  similar  observation  was 
made  with  respect  to  sum-of-products  expres¬ 
sions,  where  sum  is  MAX  and  products  are  the 
MIN  of  unary  functions  of  varying  costs  [12]. 


m.  MINIMIZATION  ALGORITHMS  FOR 
SUM  PLA’S 

We  consider  six  synthesis  methods.  The  first 
five  are  heuristic  and  are  based  on  the  direct 
cover  approach.  In  this  approach,  a  minterm  is 
first  chosen.  Next,  an  implicant  is  chosen  that 
covers  this  minterm.  The  implicant  is  then 
subtracted  from  the  function  and  the  process  is 
repeated  until  there  are  no  more  minterms. 

1.  Random  Minterm  Random  Implicant 

In  this  heuristic,  both  the  minterm  and  impli¬ 
cant  are  chosen  randomly,  with  all  choices  equally 
likely.  That  is,  among  all  minterms,  one  is 
chosen  randomly,  and  similarly  in  the  case  of 
implicants.  As  with  the  other  algorithms,  the 
chosen  implicant  is  subtracted  and  the  process 
repeated  on  the  resulting  function.  Since  no  par¬ 
ticular  characteristics  of  the  function  are  used  to 
determine  the  choice  of  minterm  or  implicant,  it 
provides  a  basis  of  comparison  for  the  next  four 
heuristics,  which  use  characteristics  of  the  func¬ 
tions  to  limit,  to  varying  degrees,  the  choices  of 
minterms  and  implicants. 

2.  Pomper  and  Armstrong 

As  in  the  Random  Minterm  Random  Impli¬ 
cant  heuristic,  this  heuristic  chooses  the  minterm 
randomly.  However,  the  implicant  is  chosen  as 


the  one,  which  when  subtracted,  drives  the  larg¬ 
est  number  of  minterms  to  0  or  don’t  care.  If 
there  is  more  than  one  such  implicant,  the  largest 
is  chosen.  If  there  are  more  than  one  largest,  the 
one  generated  first  is  chosen.  The  process  is 
repeated  until  the  function  is  completely  covered. 
A  formal  description  of  this  algorithm  is  given  in 
Appendix  I. 

Because  the  Random  Minterm  Random  Impli¬ 
cant  heuristic  can  choose  from  all  possible  covers 
of  any  function,  there  is  a  nonzero  probability 
that  it  will  find  a  minimal  realization.  However, 
if  the  minimal  realizations  are  a  small  percentage 
of  the  total  number  of  realizations,  non-minimal 
realizations  will  most  likely  result.  Our  analysis 
shows  this  clearly.  There  is  then  the  question  of 
whether  the  Pomper  and  Armstrong  heuristic  has 
a  nonzero  probability  of  finding  a  minimal  reali¬ 
zation  for  all  functions.  The  example  in  Fig.  1  is 
a  function  where  this  probability  is  0.  Here,  a 
prime  implicant  is  never  in  a  minimal  sum-of- 
products  realization,  and  so,  for  this  case,  the 
Pomper  and  Armstrong  heuristic  does  not  pro¬ 
duce  a  minimal  realization. 

3.  Besslich 

Besslich  [3]  presents  a  direct  cover  approach 
for  the  heuristic  minimization  of  multiple-valued 
logic  functions  using  minterm  weighting  and 
implicant  detecting  transformations.  In  this 
heuristic,  each  minterm  is  assigned  a  weight  that 
is  a  measure  of  the  degree  to  which  other  min¬ 
terms  cluster  around  it.  Minterms  in  the  center 
of  clusters  have  the  highest  weight  and  isolated 
minterms  the  lowest.  The  minterm  with  the 
smallest  weight  is  chosen.  Next  all  implicants 
that  cover  the  chosen  minterm  are  generated.  For 
each,  an  efficiency  factor  equal  to  the  cost  of  the 
implicant  divided  by  the  number  of  minterms  it 
covers  is  calculated  and  the  minterm  with  the 
largest  factor  is  chosen.  In  a  PLA,  all  implicants 
are  realized  with  the  same  cost  (one  column). 
Therefore,  for  this  case,  only  the  number  of  min¬ 
terms  driven  to  0  or  don’t  care  determines  which 
implicant  is  chosen.  The  basic  idea  of  the 
Besslich  heuristic  is  to  cover  the  most  isolated 
minterms  first  and  use  implicants  that  have  a  low 
cost  per  minterm  covered.  Besslich  [3]  does  not 
mention  of  how  ties  are  broken.  We  choose  to 
break  ties  among  implicants  having  the  same 
efficiency  factor  by  using  only  the  largest  impli¬ 
cants,  and  among  these  by  choosing  the  first  one 
generated.  Thus,  the  choice  of  implicants  is  made 
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in  the  same  way  as  in  the  Pomper  and 
Armstrong  heuristic.  A  formal  description  of  the 
algorithm  is  given  in  Appendix  II. 

4.  Dueck  and  Miller 

Dueck  and  Miller  [6]  present  a  heuristic  simi¬ 
lar  to  that  of  Besslich’s  with  the  intent  of 
improving  realizations  using  the  SUM  operation. 
In  this  approach,  an  isolation  factor  (IF)  is  calcu¬ 
lated  for  each  minterm  with  the  smallest  value 
(that  is,  all  1  minterms  are  considered  first;  if 
there  are  none,  then  2  minterms  are  considered, 
etc.).  The  isolation  factor  of  a  minterm  at  x  is 
inversely  proportional  to  1  plus  the  number  of 
adjacent  minterms  plus  the  number  of  logic  vari¬ 
ables  where  there  are  a  nonzero  number  of  such 
minterms.  Similar  to  the  weight  transform 
presented  by  Besslich,  the  isolation  factor  pro¬ 
vides  a  measure  of  the  degree  to  which  a  specific 
minterm  can  combine  with  other  minterms  in  the 
function  (the  correlation  is  negative,  however). 
The  minterm  with  the  highest  isolation  factor  is 
chosen.  All  implicants  that  cover  this  minterm 
are  then  generated.  For  each  of  these,  a  parame¬ 
ter  called  the  relative  break  count  (RBC)  is 
calculated.  This  provides  a  measure  of  the 
degree  to  which  the  function  is  simplified  if  the 
implicant  under  consideration  is  chosen.  The 
idea  is  to  judiciously  choose  implicants  that  make 
the  remaining  function  easy  to  realize.  Our  use 
of  "interval  literals  requires  an  adjustment  to  the 


Dueck  and  Miller  algorithm.  A  formal  descrip¬ 
tion  of  our  algorithm,  modified  to  the 
specifications  of  our  problem,  is  given  in  Appen¬ 
dix  III. 

5.  Gold 

Gold  is  a  heuristic  in  which  the  heuristic  algo¬ 
rithms  of  Pomper  and  Armstrong,  Besslich,  and 
Dueck  and  Miller  are  applied  and  the  best  realiza¬ 
tion  is  chosen.  It  was  inspired  by  the  observation 
that  these  algorithms  displayed  a  diversity  in 
realizations.  That  is,  no  single  algorithm  is  con¬ 
sistently  better  than  the  others  over  all  functions; 
there  are  classes  of  functions  where  one  algorithm 
does  better  than  the  others. 

6.  Absolute  Minimization 

An  algorithm  which  produces  the  exact 
minimal  sum-of-products  was  devised  to  compare 
with  the  results  produced  by  the  previous  four 
algorithms.  This  does  essentially  an  exhaustive 
search  over  all  possible  solutions,  starting  with 
the  fewest  number  of  product  terms.  The  algo¬ 
rithm  requires  considerable  computer  time.  A 
complete  description  of  the  algorithm  is  given  in 
Section  VI. 

7.  Summary 

Table  I  below  summarizes  the  six  algorithms. 


Algorithm 

Choice  of 

Choice  of 

Minterm 

Implicant 

1.  Random  Min.  Random  Impli. 

Random 

Random 

2.  Pomper  and  Armstrong  [13] 

Random 

Drives  most  min.  to  0 

3.  Besslich  [3] 

Smallest  weight 

Drives  most  min.  to  0 

4.  Dueck  and  Miller  [6] 

Largest  IF 

Smallest  RBC 

5.  Gold 

Best  of  2,  3,  and  4 

6.  Absolute  Minimization 

Exhaustive  search 

Table  1.  Summary  of  Minimization  Algorithms 
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IV.  COMPARISON  OF  ALGORITHM 
PERFORMANCE 


For  the  purpose  of  comparison,  we  separate 
all  four-valued,  two-variable  functions  into  17 
disjoint  classes  according  to  the  number  of  non¬ 
zero  values  in  the  function.  For  14  of  the  17 
classes,  500  random  functions  are  generated.  For 
each  function,  the  six  algorithms  are  applied  and 
the  number  of  product  terms  derived.  For  func¬ 
tions  with  less  than  three  non-zero  values,  the 
average  number  of  product  terms  can  be  calcu¬ 
lated  as  follows. 


0  non-zero  values:  There  is  only  one  (trivial) 
function,  constant  zero,  in  this  class.  This 
function  requires  no  implicants  in  its  cover, 
and  all  algorithms  produce  identical  results. 


1  non- zero  value:  There  are  48  functions  in 
this  class,  and  each  requires  uniquely  one 
implicant.  Thus,  the  average  number  of  impli¬ 
cants  required  by  all  algorithms  is  1. 


2  non-zero  values:  There  are  3^- 


such  functions.  All  require  two  implicants 
unless  there  are  two  adjacent  minterms  with 
the  same  value,  in  which  case  a  single  impli¬ 


cant  is  sufficient.  Of  the  1080  functions,  72 
can  be  covered  by  one  implicant.  Only  the 
Random  Minterm  Random  Implicant  algo¬ 
rithm  will  fail  to  find  the  minimal  cover  for 
all  functions  in  this  class.  Specifically,  func¬ 
tions  with  adjacent  identical' minterms  will  be 
incorrectly  minimized  half  the  time,  since  a 
cover  of  a  single  minterm  or  two  adjacent 
minterms  are  equally  likely.  For  a  uniformly 
distributed  set  of  functions,  the  average  and 
standard  deviation  on  the  number  of  impli¬ 
cants  produced  by  the  Random  Minterm  Ran¬ 
dom  Implicant  algorithm  are  1.97  and  0.18, 
respectively.  On  the  other  hand,  these  values 
for  all  other  algorithms  are  1.93  and  0.25, 
respectively.  Our  random  function  generation 
program  was  tested  on  this  class  and  yielded 
values  very  close  to  these,  to  within  0.18%  for 
the  average  and  to  within  3.9%  for  the  stan¬ 
dard  deviation. 

For  classes  with  a  larger  number  of  non- zero 
values,  500  functions  were  randomly  generated. 
The  functions  generated  had  no  don’t  care  values. 
However,  they  could  develop  don’t  care  values 
during  the  covering  process.  (Nota  bene:  A  value 
of  3  in  the  original  function  becomes  a  don’t  care 
when  fully  covered.)  Fig.  2  shows  the  results. 


Random  Minterm  ft  Implicant  — 


Pomper  ft  Armstrong 
Average  Number  of  Besslich 


Number  of  Nonzero  Values 

Figure  2.  The  Average  Number  of  Product  Terms  Verses  the  Number  of  Nonzero  Values 

for  the  Six  Algorithms. 
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It  can  be  seen  that  the  Random  Minterm  Ran¬ 
dom  Implicant  algorithm  does  quite  poorly. 
Choosing  the  largest  (prime)  implicant,  as  in  the 
Pomper  and  Armstrong  algorithm,  provides  a 
significant  improvement.  There  is  very  little 
difference  between  the  performance  of  the  Besslich 
and  the  Dueck  and  Miller  algorithms;  their 
curves  in  Fig.  2  are  almost  identical.  Both  algo¬ 
rithms  choose  the  most  isolated  minterm  first, 
and  this  gives  some  improvement  over  that  of  the 
Pomper  and  Armstrong  algorithm.  Gold,  taking 
advantage  of  the  small  overlap  between  the  three 
heuristic  algorithms,  provides  some  further 
improvement.  Still  further  improvement  is 
achieved  by  absolute  minimization,  shown  as  the 
lowest  curve  in  the  figure. 

Another  measure  of  comparison  of  the  algo¬ 
rithms  is  the  number  of  functions  for  which  the 
absolute  minimal  solution  is  found.  Fig.  3  shows 
how  the  six  algorithms  compare  on  this  basis. 
Here  there  is  a  much  larger  distinction  between 
the  algorithms.  As  the  number  of  non- zero 
values  increases,  the  number  of  functions  for 
which  an  optimum  cover  is  found  by  the  Random 
Minterm  Random  Implicant  algorithm  drops  off 
sharply,  to  nearly  zero.  The  Pomper  and 
Armstrong  algorithm  offers  considerable  improve¬ 


ment,  but  even  in  this  case  only  about  39%  of 
the  functions  are  completely  minimized  when  the 
number  of  non-zero  values  is  14-16.  The  Besslich 
and  Dueck  and  Miller  algorithms  are  again  very 
close  to  each  other.  Approximately,  52%  of  the 
functions  are  minimized  for  these  two  algorithms 
when  the  number  of  non-zero  values  is  14-16. 
Gold  shows  a  reasonable  improvement,  minimiz¬ 
ing  about  75%  of  the  functions  with  14-16  non¬ 
zero  values. 

Tables  II  and  HI  show  the  numerical  values 
obtained  from  the  programs.  Another  run  was 
made  using  different  seeds  for  the  random 
number  generator  streams.  The  results  from 
both  runs  were  very  similar,  indicating  that  a 
sufficiently  large  sample  set  size  was  used.  An 
examination  of  the  data  suggests  the  truth  of  the 
following. 

Conjecture;  No  more  than  10  implicants  are 
needed  in  a  minimal  sum-of-products  expres¬ 
sion  of  a  four-valued,  two-variable  function, 
where  sum  is  the  SUM  operation  and  the  pro¬ 
ducts  consist  of  the  MIN  of  a  constant  and 
interval  literals  on  the  variables. 

The  absolute  minimization  program  and  the 


Absolute  Minimum 
Gold 


Number  of  Functions  Where 
Absolute  Minimi,  is  Achieved 


Dueck  ft  Miller 
Besslich 

Pomper  ft  Armstrong 


Number  of  Nonzero  Values 

Figure  3.  The  Number  of  Functions  Where  the  Absolute  Minimal  Realization  is  Achieved 
Verses  the  Number  of  Nonzero  Values. 
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Number 

of 

Nonzero 

Values 


Random 

Minterm 

Random 

Implicant 

Pomper 

Armstrong 

Besslich 

Dueck 

Miller 

Gold 

Absolute 

Minimiza. 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.97 

1.93 

1.93 

1.93 

1.93 

1.93 

2.90 

2.79 

2.78 

2.78 

2.78 

2.78 

3.80 

3.58 

3.58 

3.57 

3.57 

3.57 

4.66 

4.35 

4.34 

4.33 

4.31 

4.31 

5.53 

5.04 

5.02 

5.00 

4.96 

4.96 

6.31 

5.66 

5.63 

5.57 

5.53 

5.52 

7.04 

6.13 

6.07 

6.03 

5,95 

5.94 

7.77 

6.67 

6.59 

6.55 

6.44 

6.40 

8.50 

7.06 

6.97 

6.91 

6.75 

6.69 

8.99 

7.34 

7.25 

7.19 

6.98 

6.92 

9.78 

7.77 

7.56 

7.57 

7.32 

7.20 

10.31 

7.90 

7.63 

7.60 

7.31 

7.17 

10.54 

7.91 

7.62 

7.60 

7.30 

7.08 

11.00 

7.81 

7.64 

7.63 

7.26 

6.97 

11.21 

7.55 

7.33 

7.24 

6.99 

6.71 

Table  11.  Average  Number  of  Product  Terms 
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random  function  generation  routines  were  written 
in  C.  The  heuristic  algorithms  were  written  in 
Pascal.  The  program  was  executed  on  a 
Hewlett-Packard  Series  9000  Model  350  worksta¬ 
tion  running  HP-UX  (HP’s  extension  of  the  UNIX 
operating  system).  The  program  took  more  than 
3  days  to  minimize  7000  randomly  generated 
functions  with  3  to  16  non-zero  values. 


V.  ABSOLUTE  MINIMIZATION. 

The  program  to  compute  the  absolute 
minimal  realization  of  a  function  is  a  search  over 
combinations  of  implicants.  As  observed  in  Sec¬ 
tion  II,  it  is  necessary  to  consider  all  implicants 
and  not  just  the  prime  implicants.  In  the  worst 
case  for  two-variable  four-valued  functions 
{f{x^,X2)  =3),  this  can  be  large,  300  implicants. 
Consider  a  function  which  has  100  implicants 
and  requires  8  product  terms  in  a  minimal  cover. 
Using  a  brute  force  enumeration  technique,  the 
proof  that  this  is  the  minimum,  requires  that  all 
combinations  of  seven  or  fewer  implicants  be 
examined.  There  are  more  than  17  billion  such 
combinations,  and  assuming  that  each  combina¬ 
tion  can  be  generated  and  examined  in  200 
microseconds  (a  value  that  is  better  than  in  our 
program),  it  would  take  more  than  944  hours,  or 
39  days,  to  examine  all  the  combinations  for  one 
function!  Clearly  this  is  not  acceptable.  However, 
using  the  lemma  given  below,  the  search  tree 
can  be  significantly  reduced. 

Definition:  U{x)  is  a  constrained  implicant  set 
of  minterm  x  for  function  /  (X),  if 
U{x)  =  /(X),  where  /(X)  is  an  impli- 

Ih)>  0 

cant  of  /  (X). 

Lemma  1:  If  U{x)  is  a  constrained  implicant  set, 
then  every  minimal  sum-of-products  expres¬ 
sion  of  /(X)  contains  at  least  one  implicant 
from  U[x). 

The  proof  follows  from  the  fact  that  at  least  one 
implicant  is  needed  from  each  of  the  k  disjoint 
sets. 

Definition:  U{x)  is  a  minimal  constrained  impli¬ 
cant  set  of  a  function  if  and  only  if 
0<  |f/(:c)|  <  |f7(y)|  for  all  yEX,  where 
t7(x)  is  a  constrained  implicant  set. 

A  minimal  constrained  implicant  set  is  useful 
in  an  efficient  search  for  the  absolute  minimal 


sum-of-products  expression  for  a  function.  The 
search  space  can  be  represented  by  a  tree  where 
each  node  corresponds  to  a  function,  and  each 
edge  to  an  implicant.  The  root  node  represents 
the  input  function,  and  all  other  nodes  represent 
functions  obtained  by  subtracting  from  the  input 
function,  the  implicants  in  the  path  (from  this 
node)  to  the  root.  If  a  function  has  ^  implicants, 
the  root  has  ^  descendents.  Further,  each  of 
these  has  at  most  (and  probably  less  than)  ^  —  1 
descendents.  However,  rather  than  consider  all 
implicants  at  the  root  node,  it  is  sufficient  to  con¬ 
sider  only  implicants  from  a  minimal  con¬ 
strained  implicant  set,  considerably  reducing  the 
search  space.  From  Lemma  1,  at  least  one  of  this 
set  must  be  in  a  minimal  sum-of-products  expres¬ 
sion.  For  example,  the  average  number  of  impli¬ 
cants  for  a  four-valued,  two-variable  function 
with  one  zero,  (obtained  from  a  randomly  gen¬ 
erated  sample  set  of  500  functions)  is  111.2, 
while  the  average  size  of  the  minimal  constrained 
implicant  set  is  only  7.9.  This  is  a  considerable 
reduction.  In  fact,  without  it  the  computation  of 
the  absolute  minimum  algorithms  would  be 
impossible. 

The  algorithm  recursively  finds  the  absolute 
minimum  cover  of  the  function  at  each  node.  Let 
G{J )  denote  the  minimum  number  of  implicants 
ne^ed  in  the  minimal  sum-of-products  expression 
of  function  /(X).  Let  tf(®)  =  {/i(X),/2(X),..., 
4(X)}  be  any  constrained  implicant  set  of  /(X). 
Let  gf{X)  =  f  {X]  —  Ii{X)  be  the  function 
obtain^  by  subtracting  /,(X)  from  /(X).  From 
the  definition  of  a  constrained  implicant  set,  it 
follows  that,  G(/ )  =  1 -f  min{G'(^i),G(5'2)v) 

ABSOLUTE  MINIMIZATION  ALGORITHM 

f  ■4—  the  input  function; 
cur_best_soln_size  <—  M  {see  below}; 
cur_soln_8ize  0; 
if  f  has  no  coverable  minterms,  then 

output  0  and  stop; 
recursively  _minimize(f); 
output  cur_best_soln_size; 
stop; 

procedure  recursively _minimize(f); 

U  some  constrained  implicant  set  of  f; 
while  ((there  exists  another  implicant  in  U) 
AND  ((cur_soln_8ize  -f-  1)  < 

(cur_best_  80ln_size)))  do  begin 

I  ■*—  the  next  implicant  in  U; 
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add  I  to  the  current  solution; 
cur_soln_size  ■«—  curjoln_size  +  1; 

if  g  has  no  coverable  minterms  then 
begin 

make  current  solution  set  the  current 
best  solution  set; 

cur_be8t_soln_size  -4—  cur  joln_size; 
end; 

else  if  ((cur_soln_size  +  1)  < 

(cur_bestjjsoln_size))  then 
begin 

recursively_minimize(g) ; 
end; 

delete  I  from  the  current  solution; 
cur_soln_size  ■*—  cur_soln_size  -  1; 
end; 

M  represents  the  starting  value  for  the  solu¬ 
tion  size  and  is  chosen  as  the  best  of  the  Pomper 
and  Armstrong,  Besslich,  and  Dueck  and  Miller 
solutions. 

The  speed  of  the  absolute  minimization  algo¬ 
rithm  depends  on  the  choice  of  the  constrained 
implicant  set.  When  U{x)  is  small,  there  are 
fewer  choices  for  an  implicant  to  cover  x,  and 
thus  less  time  is  required.  This  suggests  that  all 
possible  sets  be  scanned  and  the  smallest  chosen. 
Although  this  reduces  the  total  number  of  impli- 
cants  scanned,  the  extra  time  required  at  each 
node  may  increase  the  overall  program  execution 
time.  A  heuristic  was  used  to  decide  whether  to 
spend  this  extra  time  or  not.  In  our  implemen¬ 
tation,  U{x)  was  first  generated  by  taking  a  cov¬ 
erable  minterm  in  the  function.  An  improve¬ 
ment  was  attempted  only  if  the  number  of 
implicants  in  this  set  was  greater  than  some 
threshold  T.  Values  of  6  and  9  were  tried  with 
the  former  yielding  a  faster  program. 

VI.  CONCLUDING  REMARKS 

There  are  three  surprising  results  from  this 
study.  First,  none  of  the  three  heuristic  sum-of- 
products  minimization  algorithms  [3,6,13]  is 
clearly  superior  to  the  others.  Second,  Gold,  an 
algorithm  that  picks  the  best  of  the  three,  does 
significantly  better  than  any  particular  one.  This 
indicates  that  the  nearly  identical  average  case 
performance  is  due  to  a  diversity  of  performance 
on  individual  functions.  Third,  the  computation 
times  for  absolute  minimization  far  exce^ed  our 
pvnectatiua.  Th”®  what  b^an  as  an  investiga¬ 


tion  of  heuristic  methods  ended  as  an  investiga¬ 
tion  of  methods  to  reduce  the  search  required  by 
absolute  minimization. 
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APPENDIX  I 

POMPER  AND  ARMSTRONG  [13]  ALGORITHM 
USING  INTERVAL  LITERALS 

Orig_f  is  the  given  r-valued  n-variable  function,  f  is 
an  intermediate  function  as  it  is  covered  in  successive 
stages  of  the  algorithm.  Each  minterm  in  f  can  have  a 
value  in  the  range  [0...r]  where  r  is  the  radix.  A  value  of  - 

0  was  zero  in  orig_X  or  was  in  the  range  [l.,.(r-2)]  in  orig_f 
but  has  now  been  fully  covered, 

l,2,...,(r-l)  still  needs  to  be  covered, 
r  (don’t  care)  was  r  in  orig_f,  or  was  (r-1)  in  orig_f  but 
has  now  been  fully  covered. 

Algorithm: 

1.  f  ■* —  orig_f 

2.  If  f  has  no  minterms  in  the  range  [l...(r-l)]  then 

STOP;  else  choose  one,  randomly. 

3.  From  among  all  implicants  that  cover  Qfjniu, 
choose  the  maximal  implicant,  1. 

4.  Include  I  in  the  solution. 

5.  f  1— f  -  I 

6.  Go  to  step  2. 


Choosing  the  maximal  implicant 


An  .  implicant  can_ 

7  _ „*i_yi*2  12  .  .  . 

i  —  p  X2  . 


be 


written 


Let  the  size  of  this  implicant  be  the  number  of 
locations  where  it  is  p,  which  is  in  the  range  [l..(r- 
1)].  The  size  is, 


n  Um  -  *  m  +  1) 

l<m<n 

An  implicant  J  is  ’larger’  than  implicant  I  if  the 
size  of  J  b  larger  than  the  size  of  I.  The  maximal 
implicant  is  chosen  as  follows: 

cur_size  ■* —  -  oo 

cur_terms_covered  < OO 

for  each  implicant  I  of  f  that  completely  covers 
“min’  do  begin 

l_terms_covered  < —  the  number  of  minterms  in 
f  that  would  be  driven  to  0  or  r  (don’t  care)  if  I 
is  subtracted  from  f; 

Ljsize  < —  size  of  implicant  I; 

if  (I_terrqs_covered  >  cur_terms_covered)  OR 
((I_terms_covered  =  cur_terms_ 

covered)  AND  (I_size  >  cur_size))  then  begin 
bestAmplicant  I; 
cur_size  * —  Lsize; 

cur_terms_covered  ■* —  I_terms_covered; 
end; 
end; 


APPENDIX  II 

BESSLICH  [3]  ALGORITHM  USING  INTERVAL 
LITERALS 


Algorithm: 

Each  location  in  orig_f  corresponds  to  an  assign¬ 
ment  of  values  in  the  range  [0..(r-l)]  to  the  vari¬ 
ables  x^,X2,---,X„'  The  concatenation  of  these 
values  is  taken  as  an  index  which  represents  this 
location. 

1.  f  ■* —  orig_f 

2.  If  f  has  no  terms  in  the  range  [l...(r-l)],  then 

STOP;  else  find  the  most  isolated  minterm 
®mm 

3.  From  among  all  implicants  that  cover 
choose  the  maximal  implicant,  I. 

4.  Include  I  in  the  solution. 

5.  f  •«— f  -  I 

6.  Go  to  step  2. 


Finding  the  most  isolated  minterm 

Code  the  function  values  by  mapping  all  O’s  to  -1; 
all  r’s  to  0;  and  all  other  values  to  1.  Let  cfn  be 
the  coded  form  of  the  function  f. 
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Take  the  weight  transform  [3]  of  cfn  and  pick  the 
minterm  with  the  smallest  weight  as  follows: 

cur_lowest_wt  ■* —  oo 

for  each  term  a  in  cfn  such  that  cfn(Q!)  =  1  do 
begin 

/*  compute  the  weight  of  this  term  with  all 
other  terms  and  add  */ 

wt(a)  *-  XI  [cfn(/0)*w(a,^] 

P&fn 

if  (wt(a)  <  cur_lowe8t_wt)  then  •«—  a; 

end; 

On  exit  from  the  for  loop,  QJjjjjj,  is  chosen  as  the 
most  isolated  minterm.  In  the  loop,  w(Q!,)8)  is  com¬ 
puted  as  follows. 

Let  *„  *„_i  •••*2*1  and  '  j^j'l 

represent  the  terms  oc  and  Then,  compute 

D{a,P)=  E  Um-3m\ 

l<m<n 

and  w{a,p)  = 

Choosing  the  maximal  implicant  -  Same  as  with 
Pomper  and  Armstrong 


APPENDIX  m 

DUECK  AND  MILLER  [6]  ALGORITHM  USING 
INTERVAL  LITERALS 

Algorithm: 

1.  f  ■* —  orig_f 

2.  If  f  has  no  terms  in  the  range  [l...(r-l)],  then 
STOP;  else  find  the  most  isolated  minterm 

®inin  *** 

3.  From  among  all  implicants  that  cover  Q!n,in, 
choose  the  maximal  implicant,  I. 

4.  Include  I  in  the  solution. 

5.  f  ■< —  f  -  I,  if  f  >  I;  otherwise,  I  ■* —  r. 

6.  Go  to  step  2. 

Finding  the  most  isolated  minterm 

Let  min_val  be  the  minimum  value  of  all  terms 
in  f  that  are  in  the  range  l...(r-l)  i.e,  neither  0 
nor  don’t  care.  For  each  term  (X  in  the  function 
with  value  equal  to  min_val,  compute  the  cluster¬ 
ing  factor  (CF)  as  follows  - 

CF(ck)  =  DEA  -j-  EA^^, 


where  EAa  is  the  number  of  minterms  with 
which  oc  can  be  combined  in  an  interval  literal  and 
DEAfx  is  the  number  of  variables  (directions)  in 
which  oc  can  be  combined  with  a  nonzero  number 
of  minterms. 

Dueck  and  Miller  actually  define  the  isolation  fac¬ 
tor  (IF)  to  be  the  reciprocal  of  the  clustering  fac¬ 
tor  above  and  pick  the  minterm  with  the  highest 
IF.  The  same  result  can  be  obtained  by  picking 
the  minterm  with  the  smallest  CF,  avoiding  the 
division.  Ties  are  broken  arbitrarily.  Let  be 

the  minterm  thus  chosen. 

Choosing  the  maximal  implicant 

cur_rbc  < —  oo 

for  every  valid  implicant  I  in  f  that  completely  cov¬ 
ers  do  begin 

rbc  ■* —  calculated_rbc_of_J_in_f; 
if  (rbc  <  cur_rbc)  then  begin 
best_implicant  < —  I; 
cur_rbc  ■* —  rbc; 
end; 
end; 

Calculation  of  rbc 
rbc  < —  0; 

for  each  or  in  I  such  that  jo'l  ^  r  do  begin 
along  eacb  dimension  i  =  1  to  n  do  begin 

(if  |q:|  ^  I  l|)  OR  (if  the  neighbor  P  immedi¬ 
ately  before  oc  in  this  dimension  (if  such  P 
exists)  is  not  in  I  and  is  such  that  \P\  =  lal 
-  I)  OR  (if  the  neighbor  P  immediately  after 
OC  in  this  dimension  (if  such  P  exists)  is  not 
in  I  and  b  such  that  \P\  ==  |a|-  I)  then 
rbc  ■* —  rbc  -1; 

(if  the  neighbor  P  immediately  before  oc  in 
thb  dimension  (if  such  P  exbts)  is  not  in  I 
and  b  such  that  |y8|  =  |q!|)  OR  (if  the 
neighbor  P  immediately  after  oc  in  thb 
dimension  (if  such  P  exbts)  is  not  in  I  and  b 
such  that  |j9|=  |Qr|)  then 
rbc  ■* —  rbc  -f-  1; 

end; 

end; 
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